#We are using a .dat file (survey.dat) created from the SPSS file survey.sav  taken from SPSS Survival Manual 6th Edition Julie Pallant
#http://spss.allenandunwin.com.s3-website-ap-southeast-2.amazonaws.com/data-files.html#.Wb0vvnWP-po
#Results on a survey on well being
#We need to load the file so that we can use it in R.
sdata <- read.table("../data/survey.dat")
#Setting the column names to be that used in the dataset
colnames(sdata) <- tolower(colnames(sdata))

#Look first at partial correlation
library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
names(sdata)
##   [1] "id"         "sex"        "age"        "marital"    "child"     
##   [6] "educ"       "source"     "smoke"      "smokenum"   "op1"       
##  [11] "op2"        "op3"        "op4"        "op5"        "op6"       
##  [16] "mast1"      "mast2"      "mast3"      "mast4"      "mast5"     
##  [21] "mast6"      "mast7"      "pn1"        "pn2"        "pn3"       
##  [26] "pn4"        "pn5"        "pn6"        "pn7"        "pn8"       
##  [31] "pn9"        "pn10"       "pn11"       "pn12"       "pn13"      
##  [36] "pn14"       "pn15"       "pn16"       "pn17"       "pn18"      
##  [41] "pn19"       "pn20"       "lifsat1"    "lifsat2"    "lifsat3"   
##  [46] "lifsat4"    "lifsat5"    "pss1"       "pss2"       "pss3"      
##  [51] "pss4"       "pss5"       "pss6"       "pss7"       "pss8"      
##  [56] "pss9"       "pss10"      "sest1"      "sest2"      "sest3"     
##  [61] "sest4"      "sest5"      "sest6"      "sest7"      "sest8"     
##  [66] "sest9"      "sest10"     "m1"         "m2"         "m3"        
##  [71] "m4"         "m5"         "m6"         "m7"         "m8"        
##  [76] "m9"         "m10"        "pc1"        "pc2"        "pc3"       
##  [81] "pc4"        "pc5"        "pc6"        "pc7"        "pc8"       
##  [86] "pc9"        "pc10"       "pc11"       "pc12"       "pc13"      
##  [91] "pc14"       "pc15"       "pc16"       "pc17"       "pc18"      
##  [96] "rop2"       "rop4"       "rop6"       "rmast1"     "rmast3"    
## [101] "rmast4"     "rmast6"     "rmast7"     "rpss4"      "rpss5"     
## [106] "rpss7"      "rpss8"      "rsest3"     "rsest5"     "rsest7"    
## [111] "rsest9"     "rsest10"    "rpc1"       "rpc2"       "rpc7"      
## [116] "rpc11"      "rpc15"      "rpc16"      "toptim"     "tmast"     
## [121] "tposaff"    "tnegaff"    "tlifesat"   "tpstress"   "tslfest"   
## [126] "tmarlow"    "tpcoiss"    "agegp3"     "agegp5"     "educrec"   
## [131] "lg10negaff" "mah_1"      "coo_1"      "mah_2"
ydata<- filter(sdata,!is.na(sdata$tpcoiss), !is.na(sdata$tpstress), !is.na(sdata$toptim),!is.na(sdata$tmarlow))
head(ydata)
##    id     sex age             marital child                    educ
## 1 415 FEMALES  24  MARRIED FIRST TIME   YES COMPLETED UNDERGRADUATE
## 2   9   MALES  39 LIVING WITH PARTNER   YES COMPLETED UNDERGRADUATE
## 3 425 FEMALES  48  MARRIED FIRST TIME   YES          SOME SECONDARY
## 4 307   MALES  41           REMARRIED   YES          SOME SECONDARY
## 5 440   MALES  23              SINGLE    NO COMPLETED UNDERGRADUATE
## 6 484 FEMALES  31  MARRIED FIRST TIME   YES COMPLETED UNDERGRADUATE
##            source smoke smokenum op1 op2 op3 op4 op5 op6 mast1 mast2 mast3
## 1 LIFE IN GENERAL    NO       NA   3   2   3   2   4   2     2     4     2
## 2            WORK   YES        2   2   3   4   3   5   4     2     4     2
## 3        CHILDREN    NO       NA   3   1   3   3   3   4     3     3     2
## 4            WORK    NO        0   3   1   5   3   5   1     2     4     1
## 5            WORK    NO        0   3   2   3   2   1   3     1     4     2
## 6 LIFE IN GENERAL    NO       NA   2   2   2   2   3   4     1     3     2
##   mast4 mast5 mast6 mast7 pn1 pn2 pn3 pn4 pn5 pn6 pn7 pn8 pn9 pn10 pn11
## 1     2     4     2     3   5   5   4   4   4   5   5   5   5    5    5
## 2     3     4     2     3   4   5   4   5   3   3   3   2   5    4    4
## 3     3     3     2     2   2   1   2   2   1   1   3   1   1    1    1
## 4     1     4     1     2   5   5   3   5   2   5   5   5   5    5    5
## 5     1     2     2     2   1   1   1   1   1   1   1   2   3    1    1
## 6     2     3     1     2   2   1   1   3   1   2   1   1   1    1    3
##   pn12 pn13 pn14 pn15 pn16 pn17 pn18 pn19 pn20 lifsat1 lifsat2 lifsat3
## 1    5    5    3    5    3    5    5    4    4       5       6       5
## 2    4    4    3    2    2    2    3    4    4       7       6       5
## 3    1    1    2    2    2    1    1    2    1       7       7       7
## 4    5    4    2    5    5    5    5    3    1       7       6       7
## 5    1    1    1    1    1    1    1    1    1       3       3       4
## 6    1    1    1    1    1    1    1    1    1       2       2       2
##   lifsat4 lifsat5 pss1 pss2 pss3 pss4 pss5 pss6 pss7 pss8 pss9 pss10 sest1
## 1       4       3    3    3    4    3    4    3    3    3    3     2     4
## 2       7       5    2    2    3    5    4    3    5    4    3     3     3
## 3       6       6    1    2    2    4    4    2    4    4    2     2     4
## 4       7       6    4    3    5    5    4    3    3    4    5     3     4
## 5       3       3    2    2    3    2    3    2    4    3    3     3     2
## 6       2       2    1    1    3    4    3    2    4    3    2     2     2
##   sest2 sest3 sest4 sest5 sest6 sest7 sest8 sest9 sest10 m1 m2 m3 m4 m5 m6
## 1     4     1     4     1     4     2     4     3      3  1  0  1  1  0  0
## 2     4     1     4     1     4     3     3     1      3  1  1  1  0  1  0
## 3     3     3     3     1     3     2     3     2      2  1  1  1  1  1  0
## 4     4     1     4     1     4     1     4     1      1  1  1  0  0  0  0
## 5     2     2     2     4     2     3     2     3      2  1  0  1  1  0  0
## 6     3     2     3     2     2     3     2     3      3  1  1  1  1  1  1
##   m7 m8 m9 m10 pc1 pc2 pc3 pc4 pc5 pc6 pc7 pc8 pc9 pc10 pc11 pc12 pc13
## 1  1  0  0   0   4   4   3   3   4   4   4   2   3    2    4    4    2
## 2  0  1  0   0   4   3   3   4   3   2   5   1   1    1    5    2    4
## 3  1  1  0   1   3   3   4   2   3   3   3   3   3    3    4    3    2
## 4  0  1  0   0   2   2   4   5   4   4   4   3   4    4    1    3    2
## 5  0  1  0   0   1   2   2   1   2   1   2   2   3    1    3    2    1
## 6  1  1  1   1   1   3   4   3   2   1   4   2   3    3    3    1    4
##   pc14 pc15 pc16 pc17 pc18 rop2 rop4 rop6 rmast1 rmast3 rmast4 rmast6
## 1    2    3    3    4    4    4    4    4      3      3      3      3
## 2    2    4    2    2    2    3    3    2      3      3      2      3
## 3    2    4    4    2    2    5    3    2      2      3      2      3
## 4    2    3    2    4    2    5    3    5      3      4      4      4
## 5    3    2    1    1    2    4    4    3      4      3      4      3
## 6    1    4    4    1    1    4    4    2      4      3      3      4
##   rmast7 rpss4 rpss5 rpss7 rpss8 rsest3 rsest5 rsest7 rsest9 rsest10 rpc1
## 1      2     3     2     3     3      4      4      3      2       2    2
## 2      2     1     2     1     2      4      4      2      4       2    2
## 3      3     2     2     2     2      2      4      3      3       3    3
## 4      3     1     2     3     2      4      4      4      4       4    4
## 5      3     4     3     2     3      3      1      2      2       3    5
## 6      3     2     3     2     3      3      3      2      2       2    5
##   rpc2 rpc7 rpc11 rpc15 rpc16 toptim tmast tposaff tnegaff tlifesat
## 1    2    2     2     3     3     22    22      49      39       23
## 2    3    1     1     2     4     19    21      35      35       30
## 3    3    3     2     2     2     19    19      15      14       33
## 4    4    2     5     3     4     26    26      49      36       33
## 5    4    4     3     4     5     18    23      12      11       16
## 6    3    2     3     2     2     17    23      14      12       10
##   tpstress tslfest tmarlow tpcoiss  agegp3  agegp5
## 1       29      35       4      51 18 - 29 18 - 24
## 2       22      34       5      40 30 - 44 33 - 40
## 3       19      31       8      47     45+ 41 - 49
## 4       31      40       3      63 30 - 44 41 - 49
## 5       27      21       4      46 18 - 29 18 - 24
## 6       21      24      10      43 30 - 44 25 - 32
##                        educrec lg10negaff     mah_1        coo_1    mah_2
## 1      completed undergrad uni   1.623249 0.9580162 4.534498e-05 18.10064
## 2      completed undergrad uni   1.544068 3.6407533 1.585239e-02 14.48201
## 3 did not complete high school   1.146128 1.3074179 1.404248e-02 14.21429
## 4 did not complete high school   1.556303 1.3141761 5.403849e-03 14.16941
## 5      completed undergrad uni   1.041393 2.7188874 3.755638e-04 13.70297
## 6      completed undergrad uni   1.079181 3.7542989 1.351416e-02 13.56551
library(ppcor)
## Loading required package: MASS
## 
## Attaching package: 'MASS'
## The following object is masked from 'package:dplyr':
## 
##     select
#perception of control and stress controlling for social desirability
spcor.test(ydata$tpcoiss, ydata$tpstress, ydata$tmarlow)
##     estimate      p.value statistic   n gp  Method
## 1 -0.5269534 1.862222e-31 -12.69154 422  1 pearson
#perception of control and optimism controlling for social desirability
spcor.test(ydata$tpcoiss, ydata$toptim, ydata$tmarlow)
##    estimate      p.value statistic   n gp  Method
## 1 0.4830465 5.348979e-26  11.29257 422  1 pearson
#stress and optimism controlling for social desirability
spcor.test(ydata$tpstress, ydata$toptim, ydata$tmarlow)
##     estimate     p.value statistic   n gp  Method
## 1 -0.4412444 1.73877e-21 -10.06483 422  1 pearson
#Look at differences in scores for gender
#Conduct Levene's test for homogeneity of variance in library car
library(car)
## Loading required package: carData
## 
## Attaching package: 'car'
## The following object is masked from 'package:dplyr':
## 
##     recode
leveneTest(tpstress ~ sex, data=ydata)
## Levene's Test for Homogeneity of Variance (center = median)
##        Df F value Pr(>F)
## group   1   1.678 0.1959
##       420
#Conduct the t-test
#You can use the var.equal = TRUE option to specify equal variances and a pooled variance estimate
t.test(tpstress~sex,var.equal=FALSE,data=ydata)#Variances are not equal
## 
##  Welch Two Sample t-test
## 
## data:  tpstress by sex
## t = 3.0956, df = 402.98, p-value = 0.002101
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  0.635108 2.845362
## sample estimates:
## mean in group FEMALES   mean in group MALES 
##              27.46091              25.72067
#Conduct Levene's test for homogeneity of variance in library car
leveneTest(tpcoiss ~ sex, data=ydata)
## Levene's Test for Homogeneity of Variance (center = median)
##        Df F value Pr(>F)
## group   1  2.2499 0.1344
##       420
#Conduct the t-test
#You can use the var.equal = TRUE option to specify equal variances and a pooled variance estimate
t.test(tpcoiss~sex,var.equal=FALSE,data=ydata)
## 
##  Welch Two Sample t-test
## 
## data:  tpcoiss by sex
## t = -2.1841, df = 401.99, p-value = 0.02953
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -4.8177980 -0.2533563
## sample estimates:
## mean in group FEMALES   mean in group MALES 
##              59.49794              62.03352
#Conduct Levene's test for homogeneity of variance in library car
leveneTest(toptim ~ sex, data=ydata)
## Levene's Test for Homogeneity of Variance (center = median)
##        Df F value  Pr(>F)  
## group   1  4.5018 0.03444 *
##       420                  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#Conduct the t-test
#You can use the var.equal = TRUE option to specify equal variances and a pooled variance estimate
t.test(toptim~sex,var.equal=TRUE,data=ydata)
## 
##  Two Sample t-test
## 
## data:  toptim by sex
## t = 0.40731, df = 420, p-value = 0.684
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -0.6780665  1.0325278
## sample estimates:
## mean in group FEMALES   mean in group MALES 
##              22.23868              22.06145
#Conduct Levene's test for homogeneity of variance in library car
leveneTest(tmarlow ~ sex, data=sdata)
## Levene's Test for Homogeneity of Variance (center = median)
##        Df F value Pr(>F)
## group   1  0.0132 0.9087
##       431
#Conduct the t-test
#You can use the var.equal = TRUE option to specify equal variances and a pooled variance estimate
t.test(tmarlow~sex,var.equal=FALSE,data=ydata)
## 
##  Welch Two Sample t-test
## 
## data:  tmarlow by sex
## t = 1.9176, df = 388.54, p-value = 0.05589
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -0.009565282  0.766720947
## sample estimates:
## mean in group FEMALES   mean in group MALES 
##              5.456790              5.078212
#install the library userfriendlyscience and load it using the library command
#it has a really nice one-way anova function that provides
#nice summary output

library(userfriendlyscience) 

#run a one-way anova test using the correct post-hoc test Tukey in our case
#Use Games-Howell for unequal variances
one.way <- oneway(ydata$agegp3, y = ydata$toptim, posthoc = 'Tukey') 
#printout a summary of the anova 
one.way 
## ### Oneway Anova for y=toptim and x=agegp3 (groups: 18 - 29, 30 - 44, 45+)
## 
## Omega squared: 95% CI = [0; .06], point estimate = .02
## Eta Squared: 95% CI = [0; .05], point estimate = .02
## 
##                                      SS  Df    MS    F    p
## Between groups (error + effect)  196.06   2 98.03 5.13 .006
## Within groups (error only)      8003.66 419  19.1          
## 
## 
## ### Post hoc test: Tukey
## 
##                 diff lwr   upr  p adj
## 30 - 44-18 - 29 0.66 -0.54 1.86 .403 
## 45+-18 - 29     1.69 0.44  2.94 .004 
## 45+-30 - 44     1.04 -0.2  2.27 .121
one.way <- oneway(ydata$agegp3, y = ydata$tpstress, posthoc = 'Tukey') 
#printout a summary of the anova 
one.way
## ### Oneway Anova for y=tpstress and x=agegp3 (groups: 18 - 29, 30 - 44, 45+)
## 
## Omega squared: 95% CI = [NA; .04], point estimate = .01
## Eta Squared: 95% CI = [0; .03], point estimate = .01
## 
##                                       SS  Df    MS    F    p
## Between groups (error + effect)   160.54   2 80.27 2.36 .096
## Within groups (error only)      14262.02 419 34.04          
## 
## 
## ### Post hoc test: Tukey
## 
##                 diff  lwr   upr  p adj
## 30 - 44-18 - 29 -1.01 -2.61 0.59 .301 
## 45+-18 - 29     -1.5  -3.16 0.17 .089 
## 45+-30 - 44     -0.49 -2.14 1.16 .767
one.way <- oneway(ydata$agegp3, y = ydata$tpcoiss, posthoc = 'Tukey') 
#printout a summary of the anova 
one.way
## ### Oneway Anova for y=tpcoiss and x=agegp3 (groups: 18 - 29, 30 - 44, 45+)
## 
## Omega squared: 95% CI = [.01; .08], point estimate = .04
## Eta Squared: 95% CI = [.01; .08], point estimate = .04
## 
##                                       SS  Df      MS    F     p
## Between groups (error + effect)  2626.31   2 1313.16 9.47 <.001
## Within groups (error only)      58102.91 419  138.67           
## 
## 
## ### Post hoc test: Tukey
## 
##                 diff lwr   upr  p adj
## 30 - 44-18 - 29 1.34 -1.89 4.57 .592 
## 45+-18 - 29     5.97 2.61  9.34 <.001
## 45+-30 - 44     4.63 1.3   7.96 .003
one.way <- oneway(ydata$agegp3, y = ydata$tmarlow, posthoc = 'Tukey') 
#printout a summary of the anova 
one.way
## ### Oneway Anova for y=tmarlow and x=agegp3 (groups: 18 - 29, 30 - 44, 45+)
## 
## Omega squared: 95% CI = [.02; .1], point estimate = .05
## Eta Squared: 95% CI = [.02; .09], point estimate = .06
## 
##                                      SS  Df    MS     F     p
## Between groups (error + effect)   95.21   2 47.61 12.32 <.001
## Within groups (error only)      1618.76 419  3.86            
## 
## 
## ### Post hoc test: Tukey
## 
##                 diff lwr  upr  p adj
## 30 - 44-18 - 29 0.79 0.25 1.33 .002 
## 45+-18 - 29     1.15 0.58 1.71 <.001
## 45+-30 - 44     0.35 -0.2 0.91 .297

Build the linear regression models

#Baseline model optimism and social desirability as predictors
model1=lm(ydata$tpcoiss~ydata$toptim+ydata$tmarlow)
summary(model1)
## 
## Call:
## lm(formula = ydata$tpcoiss ~ ydata$toptim + ydata$tmarlow)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -27.5175  -6.4859   0.4332   6.4263  26.2576 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    23.7960     2.6606   8.944  < 2e-16 ***
## ydata$toptim    1.3249     0.1104  11.998  < 2e-16 ***
## ydata$tmarlow   1.3996     0.2415   5.795 1.35e-08 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 9.921 on 419 degrees of freedom
## Multiple R-squared:  0.3209, Adjusted R-squared:  0.3176 
## F-statistic: 98.98 on 2 and 419 DF,  p-value: < 2.2e-16
plot(model1)

#Check assumptions
#Cooks distance
cooks.distance(model1)
##            1            2            3            4            5 
## 6.506115e-04 3.129057e-03 5.087877e-03 8.452198e-06 9.240629e-04 
##            6            7            8            9           10 
## 2.161413e-02 2.048926e-02 1.193240e-02 1.939681e-02 3.356655e-02 
##           11           12           13           14           15 
## 2.146690e-02 1.566864e-02 1.287334e-09 1.277584e-02 2.627391e-03 
##           16           17           18           19           20 
## 4.739549e-03 5.190616e-03 9.023535e-04 5.028507e-03 9.223243e-04 
##           21           22           23           24           25 
## 5.089454e-04 8.069809e-03 9.064434e-04 6.878265e-04 2.129152e-03 
##           26           27           28           29           30 
## 6.964412e-04 4.149586e-03 1.524208e-04 5.381614e-04 1.224734e-03 
##           31           32           33           34           35 
## 5.821577e-02 2.108022e-04 4.576823e-05 7.420428e-05 1.952397e-03 
##           36           37           38           39           40 
## 6.069657e-04 1.150524e-03 9.023535e-04 2.875893e-04 8.964110e-03 
##           41           42           43           44           45 
## 2.513665e-03 4.173919e-03 1.593331e-03 3.366332e-04 2.249850e-03 
##           46           47           48           49           50 
## 2.567601e-04 4.621626e-06 5.983472e-05 2.968610e-04 4.357802e-03 
##           51           52           53           54           55 
## 7.879206e-03 8.223513e-03 2.919133e-04 1.551056e-04 2.820941e-03 
##           56           57           58           59           60 
## 2.721381e-05 7.439357e-03 1.214220e-03 2.455651e-04 1.955042e-02 
##           61           62           63           64           65 
## 8.217066e-03 2.531394e-02 2.938392e-05 3.485609e-03 5.595200e-04 
##           66           67           68           69           70 
## 3.067932e-03 1.221729e-04 3.111358e-04 4.399424e-04 1.178932e-03 
##           71           72           73           74           75 
## 1.153362e-03 6.837250e-06 6.168852e-03 4.165562e-03 6.094952e-04 
##           76           77           78           79           80 
## 1.845418e-02 2.968610e-04 2.848670e-04 6.064830e-04 2.401268e-03 
##           81           82           83           84           85 
## 9.691870e-04 6.827054e-04 1.021566e-03 1.001439e-04 3.270077e-04 
##           86           87           88           89           90 
## 1.723499e-03 9.107545e-07 6.971367e-03 9.013481e-04 1.658224e-04 
##           91           92           93           94           95 
## 1.402775e-03 4.077905e-03 6.812020e-04 7.180430e-03 8.656207e-05 
##           96           97           98           99          100 
## 6.582227e-05 4.193576e-05 3.458286e-03 3.127371e-05 2.264226e-04 
##          101          102          103          104          105 
## 1.392142e-03 1.354924e-03 2.522240e-05 2.631634e-03 1.063894e-06 
##          106          107          108          109          110 
## 1.924593e-03 1.060526e-02 8.806716e-05 3.650248e-03 5.889678e-05 
##          111          112          113          114          115 
## 3.010354e-05 5.381614e-04 3.067745e-04 8.517758e-03 2.963168e-03 
##          116          117          118          119          120 
## 1.176251e-03 2.936065e-04 1.650252e-03 1.656681e-03 3.342240e-05 
##          121          122          123          124          125 
## 2.066646e-03 4.975720e-03 6.012884e-03 1.763650e-05 1.274406e-04 
##          126          127          128          129          130 
## 1.487809e-04 4.591151e-05 2.235437e-03 1.025210e-04 1.402666e-03 
##          131          132          133          134          135 
## 1.513572e-03 2.295381e-04 4.692662e-03 5.194474e-03 1.632979e-03 
##          136          137          138          139          140 
## 1.624723e-04 1.198991e-04 3.583611e-03 4.056554e-04 3.335224e-04 
##          141          142          143          144          145 
## 3.797979e-05 4.819077e-03 2.773353e-03 1.499621e-04 5.689831e-05 
##          146          147          148          149          150 
## 2.862607e-03 9.999551e-03 4.909010e-03 7.464769e-04 9.777531e-05 
##          151          152          153          154          155 
## 4.087176e-03 3.036948e-05 1.757349e-04 7.157722e-04 4.139358e-04 
##          156          157          158          159          160 
## 3.495057e-03 1.924870e-06 1.350806e-06 3.964415e-04 1.039854e-06 
##          161          162          163          164          165 
## 1.091038e-03 1.372123e-04 3.445105e-04 2.179824e-03 1.175263e-03 
##          166          167          168          169          170 
## 6.338430e-03 9.987183e-05 7.884016e-04 2.073466e-03 1.554853e-04 
##          171          172          173          174          175 
## 1.198991e-04 4.994177e-03 1.380328e-03 4.696045e-03 2.464019e-03 
##          176          177          178          179          180 
## 5.179151e-03 3.382006e-03 1.468342e-07 6.529841e-03 1.144402e-02 
##          181          182          183          184          185 
## 5.708992e-03 6.957986e-03 8.605417e-05 4.538805e-04 1.210417e-03 
##          186          187          188          189          190 
## 1.170647e-03 1.397110e-03 2.642224e-03 3.434621e-06 1.053473e-04 
##          191          192          193          194          195 
## 6.040439e-03 1.021492e-05 2.621312e-05 3.058127e-04 1.106836e-06 
##          196          197          198          199          200 
## 3.084810e-04 3.681709e-04 7.158949e-06 7.232093e-04 5.107740e-05 
##          201          202          203          204          205 
## 1.677699e-04 5.885062e-04 2.053106e-04 7.578335e-04 4.736659e-04 
##          206          207          208          209          210 
## 4.538805e-04 1.919009e-02 5.959021e-06 1.174626e-03 3.559009e-04 
##          211          212          213          214          215 
## 1.179592e-07 2.617542e-03 1.439938e-03 5.375908e-04 6.156527e-03 
##          216          217          218          219          220 
## 1.949562e-04 2.161649e-04 3.375625e-04 2.521972e-03 3.342415e-04 
##          221          222          223          224          225 
## 9.711773e-05 1.052173e-03 2.070093e-03 1.642951e-03 1.752558e-04 
##          226          227          228          229          230 
## 2.359908e-04 1.380037e-05 5.696251e-04 2.503399e-05 3.978605e-03 
##          231          232          233          234          235 
## 3.412849e-03 7.974683e-03 7.578335e-04 4.340043e-06 4.957295e-05 
##          236          237          238          239          240 
## 1.072370e-04 1.524910e-03 3.167750e-08 7.816887e-04 3.989980e-04 
##          241          242          243          244          245 
## 4.328349e-03 1.195432e-03 2.627620e-02 5.271424e-03 5.104507e-03 
##          246          247          248          249          250 
## 2.922495e-03 9.186979e-03 4.017850e-04 1.271078e-03 6.135397e-07 
##          251          252          253          254          255 
## 8.145769e-04 1.688660e-03 7.337091e-05 3.999015e-05 1.078057e-04 
##          256          257          258          259          260 
## 6.388485e-04 3.409956e-04 4.957295e-05 1.792341e-04 3.701106e-04 
##          261          262          263          264          265 
## 8.368294e-05 2.935490e-03 7.674842e-03 1.131683e-03 1.280714e-04 
##          266          267          268          269          270 
## 5.598536e-05 2.223945e-03 4.416331e-03 1.857985e-03 1.672161e-02 
##          271          272          273          274          275 
## 6.659797e-03 2.483211e-05 1.966412e-03 1.061073e-04 1.129694e-03 
##          276          277          278          279          280 
## 8.682361e-03 6.418422e-05 7.193885e-03 3.001158e-03 3.342110e-04 
##          281          282          283          284          285 
## 1.372123e-04 4.290792e-03 8.865922e-06 1.651619e-03 2.033704e-04 
##          286          287          288          289          290 
## 5.499787e-04 7.280757e-04 1.439891e-03 2.904122e-04 1.256372e-03 
##          291          292          293          294          295 
## 1.294287e-03 2.284546e-03 2.700759e-02 1.490000e-03 4.149014e-04 
##          296          297          298          299          300 
## 3.109561e-03 4.616497e-04 9.050319e-06 4.594888e-07 7.625216e-03 
##          301          302          303          304          305 
## 2.566532e-03 1.649468e-03 6.496286e-06 1.914557e-05 2.785114e-03 
##          306          307          308          309          310 
## 4.364782e-04 1.757705e-03 2.343748e-03 2.867822e-03 4.054453e-06 
##          311          312          313          314          315 
## 9.264442e-04 1.746174e-03 1.452083e-03 9.805261e-04 1.584094e-05 
##          316          317          318          319          320 
## 1.008374e-04 1.256380e-05 3.010354e-05 4.532669e-04 5.594561e-04 
##          321          322          323          324          325 
## 1.823669e-03 1.941765e-05 2.936065e-04 8.714722e-05 6.589561e-04 
##          326          327          328          329          330 
## 7.834848e-03 2.967029e-05 1.210417e-03 3.524980e-03 1.711533e-03 
##          331          332          333          334          335 
## 7.967586e-04 9.576450e-04 6.905434e-05 6.421580e-05 1.571377e-04 
##          336          337          338          339          340 
## 4.793292e-05 8.724819e-04 3.883082e-04 1.658224e-04 7.986904e-03 
##          341          342          343          344          345 
## 6.088274e-04 8.053373e-04 2.427713e-03 2.391558e-03 2.106339e-03 
##          346          347          348          349          350 
## 1.716402e-04 2.674903e-03 1.605706e-04 2.612864e-03 1.126339e-04 
##          351          352          353          354          355 
## 3.366332e-04 1.137714e-04 1.818245e-04 6.503035e-04 9.047490e-03 
##          356          357          358          359          360 
## 7.790803e-04 8.288269e-04 9.748803e-05 6.492991e-07 1.508200e-03 
##          361          362          363          364          365 
## 5.473312e-04 6.387701e-03 3.991669e-03 1.119548e-03 1.642899e-03 
##          366          367          368          369          370 
## 1.519985e-02 1.620073e-05 1.062820e-03 1.780379e-04 2.819749e-03 
##          371          372          373          374          375 
## 1.501189e-03 3.294622e-04 9.732599e-06 3.213351e-03 3.193184e-03 
##          376          377          378          379          380 
## 3.661866e-04 1.492295e-03 6.734092e-04 1.018918e-06 1.862339e-04 
##          381          382          383          384          385 
## 3.255562e-03 8.866627e-06 4.400965e-03 8.037937e-04 5.023406e-05 
##          386          387          388          389          390 
## 1.319660e-02 1.553321e-03 4.840320e-04 1.572066e-03 6.204326e-05 
##          391          392          393          394          395 
## 1.231578e-04 7.545399e-05 1.878922e-03 1.025210e-04 1.466992e-03 
##          396          397          398          399          400 
## 3.367960e-06 5.141552e-05 2.602819e-03 2.461654e-03 3.410122e-05 
##          401          402          403          404          405 
## 5.132774e-04 4.397996e-05 7.133846e-05 8.613094e-03 2.990435e-03 
##          406          407          408          409          410 
## 5.617803e-04 2.694865e-07 1.408026e-03 1.080433e-02 1.438883e-03 
##          411          412          413          414          415 
## 4.256345e-05 1.779525e-02 1.242404e-03 3.588363e-04 1.279044e-06 
##          416          417          418          419          420 
## 1.117406e-04 3.575873e-04 2.949043e-04 4.768409e-04 7.689343e-04 
##          421          422 
## 4.591151e-05 4.500562e-05
#Create histogram
 
#A density plot of the residuals
plot(density(resid(model1))) 

#Create a QQ plotqqPlot(model, main="QQ Plot") #qq plot for studentized resid 
leveragePlots(model1) # leverage plots

#Collinearity
vif(model1)
##  ydata$toptim ydata$tmarlow 
##      1.015801      1.015801
sqrt(vif(model1)) 
##  ydata$toptim ydata$tmarlow 
##       1.00787       1.00787
#Second model adding in stress
model2=lm(ydata$tpcoiss~ydata$toptim+ydata$tmarlow+ydata$tpstress)
summary(model2)
## 
## Call:
## lm(formula = ydata$tpcoiss ~ ydata$toptim + ydata$tmarlow + ydata$tpstress)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -34.253  -5.034   0.271   5.701  28.070 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    58.55471    4.40587  13.290  < 2e-16 ***
## ydata$toptim    0.84267    0.11264   7.481 4.38e-13 ***
## ydata$tmarlow   0.98044    0.22400   4.377 1.52e-05 ***
## ydata$tpstress -0.81767    0.08663  -9.438  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 9.019 on 418 degrees of freedom
## Multiple R-squared:  0.4402, Adjusted R-squared:  0.4361 
## F-statistic: 109.5 on 3 and 418 DF,  p-value: < 2.2e-16
plot(model2)

#Check assumptions
#Cooks distance
cooks.distance(model2)
##            1            2            3            4            5 
## 4.386133e-04 1.122128e-02 1.841839e-02 8.104241e-04 1.677983e-03 
##            6            7            8            9           10 
## 3.951943e-02 1.278646e-02 1.858489e-02 1.124209e-02 1.319163e-02 
##           11           12           13           14           15 
## 2.095067e-02 3.536274e-02 2.078877e-04 1.510310e-02 4.006136e-04 
##           16           17           18           19           20 
## 2.818643e-03 5.081900e-03 2.784382e-02 6.187649e-04 9.385656e-04 
##           21           22           23           24           25 
## 1.692352e-04 2.913993e-02 2.425825e-03 1.963788e-03 1.334735e-03 
##           26           27           28           29           30 
## 2.386220e-05 5.446850e-03 1.513688e-03 9.679650e-07 3.836379e-04 
##           31           32           33           34           35 
## 3.596554e-02 6.542349e-04 1.963583e-04 4.635160e-04 1.006445e-03 
##           36           37           38           39           40 
## 3.779544e-04 6.207859e-04 2.893225e-03 1.518792e-03 6.252466e-03 
##           41           42           43           44           45 
## 4.803713e-04 3.623659e-03 1.038313e-02 2.437726e-03 5.997050e-04 
##           46           47           48           49           50 
## 6.230510e-03 1.268997e-04 1.895822e-03 1.023493e-03 5.358517e-03 
##           51           52           53           54           55 
## 6.697990e-02 4.617651e-02 7.650889e-04 4.571700e-03 3.066488e-03 
##           56           57           58           59           60 
## 3.898795e-04 6.702801e-03 9.716547e-04 8.870216e-05 7.531584e-03 
##           61           62           63           64           65 
## 1.070167e-02 1.516978e-02 5.397230e-05 4.659988e-04 8.643770e-03 
##           66           67           68           69           70 
## 1.326955e-03 4.032575e-03 7.038843e-04 6.651905e-04 3.143702e-04 
##           71           72           73           74           75 
## 4.152421e-04 1.621811e-05 7.732975e-04 1.689876e-02 1.567265e-05 
##           76           77           78           79           80 
## 1.402302e-02 6.050244e-04 3.492398e-04 5.618475e-04 1.600254e-03 
##           81           82           83           84           85 
## 3.831070e-06 3.579970e-04 3.061522e-04 1.108126e-03 1.156157e-03 
##           86           87           88           89           90 
## 1.376175e-03 3.209951e-04 3.646189e-03 2.361536e-04 3.462976e-04 
##           91           92           93           94           95 
## 8.202307e-05 1.680604e-04 7.306675e-04 4.747362e-03 7.856356e-05 
##           96           97           98           99          100 
## 1.008669e-05 6.809411e-04 1.238707e-02 2.383199e-04 6.062445e-04 
##          101          102          103          104          105 
## 1.233405e-03 8.383219e-04 2.944082e-04 4.329760e-04 1.011100e-03 
##          106          107          108          109          110 
## 1.938113e-03 7.450220e-03 1.025415e-04 2.037185e-03 1.223689e-04 
##          111          112          113          114          115 
## 5.507892e-06 3.817765e-04 3.338964e-05 7.695163e-03 2.624638e-03 
##          116          117          118          119          120 
## 9.518047e-04 6.441381e-05 1.690502e-03 1.098883e-03 6.995576e-04 
##          121          122          123          124          125 
## 5.990789e-03 1.280055e-02 3.799758e-03 6.759402e-05 6.630600e-05 
##          126          127          128          129          130 
## 2.479354e-04 3.051437e-05 5.508596e-04 1.585305e-04 1.050701e-03 
##          131          132          133          134          135 
## 7.095943e-04 8.935817e-04 6.350803e-04 4.076625e-02 1.429977e-02 
##          136          137          138          139          140 
## 1.206158e-03 5.763347e-04 2.135666e-03 8.705416e-09 1.884086e-03 
##          141          142          143          144          145 
## 4.320267e-05 3.078351e-03 2.194644e-03 1.164725e-05 6.621026e-05 
##          146          147          148          149          150 
## 1.402424e-02 6.955861e-03 3.331298e-03 1.002329e-03 1.530437e-04 
##          151          152          153          154          155 
## 2.578946e-03 7.411840e-06 4.302432e-04 3.132903e-03 5.284563e-05 
##          156          157          158          159          160 
## 1.908219e-03 1.067378e-03 5.755838e-04 6.092257e-04 5.196341e-05 
##          161          162          163          164          165 
## 8.455585e-04 3.543527e-04 1.318425e-04 1.881276e-03 9.157110e-04 
##          166          167          168          169          170 
## 6.946290e-03 5.300279e-05 2.127500e-05 1.295259e-03 7.237538e-05 
##          171          172          173          174          175 
## 8.855640e-05 4.740295e-03 1.064043e-03 3.319243e-03 3.817172e-04 
##          176          177          178          179          180 
## 4.436219e-03 2.669174e-03 9.347653e-04 6.462383e-03 1.053434e-02 
##          181          182          183          184          185 
## 5.679416e-03 4.930282e-03 4.027091e-04 4.504617e-04 1.643329e-03 
##          186          187          188          189          190 
## 4.818361e-05 1.592021e-03 2.571586e-03 4.890655e-05 2.127256e-03 
##          191          192          193          194          195 
## 6.103995e-03 1.201321e-04 1.864803e-04 6.200510e-05 1.986242e-04 
##          196          197          198          199          200 
## 2.126468e-04 1.579025e-05 2.576602e-05 2.475584e-04 1.012185e-04 
##          201          202          203          204          205 
## 8.324393e-05 1.107312e-03 3.054095e-07 2.217814e-04 4.681305e-04 
##          206          207          208          209          210 
## 1.115400e-03 1.223851e-02 5.757810e-06 4.472341e-04 8.394730e-04 
##          211          212          213          214          215 
## 9.312336e-06 3.066385e-03 4.856027e-04 5.113950e-04 4.688919e-03 
##          216          217          218          219          220 
## 8.904989e-04 2.971711e-04 1.265362e-04 2.381366e-03 4.750834e-06 
##          221          222          223          224          225 
## 2.481636e-04 1.304818e-04 1.020656e-03 1.199490e-02 2.910689e-04 
##          226          227          228          229          230 
## 1.204419e-04 7.463380e-06 7.239221e-04 1.370971e-05 3.586538e-03 
##          231          232          233          234          235 
## 1.352064e-03 7.742845e-03 3.192741e-05 9.414624e-06 1.400134e-05 
##          236          237          238          239          240 
## 2.247768e-04 2.751394e-03 4.638816e-06 1.092434e-03 4.336220e-05 
##          241          242          243          244          245 
## 5.012996e-03 1.846043e-03 2.276583e-02 3.502830e-03 3.373733e-03 
##          246          247          248          249          250 
## 1.637806e-03 7.437340e-03 6.446441e-05 1.249926e-03 1.628093e-05 
##          251          252          253          254          255 
## 2.089036e-03 8.269490e-04 6.272452e-04 1.555396e-04 5.870607e-05 
##          256          257          258          259          260 
## 1.277570e-03 8.920941e-05 1.653760e-03 5.820557e-06 2.716542e-04 
##          261          262          263          264          265 
## 2.781698e-05 2.372743e-03 7.082099e-03 9.568469e-04 5.097000e-06 
##          266          267          268          269          270 
## 2.253401e-04 1.783613e-03 3.420871e-03 1.023987e-03 9.543107e-03 
##          271          272          273          274          275 
## 8.586335e-03 1.260799e-05 4.736527e-03 1.722811e-07 5.328438e-05 
##          276          277          278          279          280 
## 9.064871e-03 2.146134e-05 5.317210e-03 1.760873e-03 1.905966e-05 
##          281          282          283          284          285 
## 4.640265e-05 4.380822e-03 1.093820e-04 9.877358e-04 1.344988e-04 
##          286          287          288          289          290 
## 1.585704e-04 3.916155e-04 7.051611e-04 1.610758e-04 4.548479e-04 
##          291          292          293          294          295 
## 4.538674e-04 1.201836e-03 2.143479e-02 1.205541e-03 3.031590e-04 
##          296          297          298          299          300 
## 8.962552e-04 1.114583e-04 1.393493e-04 5.159693e-06 6.983783e-03 
##          301          302          303          304          305 
## 2.251318e-03 2.902743e-03 2.288025e-04 6.227432e-04 2.346624e-03 
##          306          307          308          309          310 
## 4.089628e-05 1.387249e-03 2.760761e-03 2.495428e-03 3.530648e-05 
##          311          312          313          314          315 
## 3.380011e-04 5.835935e-04 2.496055e-03 2.774031e-04 2.308443e-04 
##          316          317          318          319          320 
## 1.686315e-04 1.143373e-03 9.677548e-05 1.683532e-03 7.132650e-04 
##          321          322          323          324          325 
## 1.353526e-03 5.250248e-06 1.256242e-04 7.959785e-05 4.593538e-04 
##          326          327          328          329          330 
## 5.517579e-03 1.932397e-05 5.013628e-04 3.298492e-03 6.990907e-04 
##          331          332          333          334          335 
## 4.655428e-05 3.522176e-04 2.364236e-05 2.207084e-05 3.881046e-05 
##          336          337          338          339          340 
## 1.627580e-05 1.514174e-03 1.055266e-04 1.433574e-04 6.123495e-03 
##          341          342          343          344          345 
## 2.574264e-04 3.304918e-04 6.589807e-04 8.551661e-04 1.323455e-03 
##          346          347          348          349          350 
## 5.606057e-06 2.238127e-03 5.255332e-05 1.280023e-03 1.999500e-05 
##          351          352          353          354          355 
## 5.766521e-05 6.914136e-05 6.995432e-05 6.633236e-04 2.361188e-02 
##          356          357          358          359          360 
## 4.722831e-04 5.031333e-03 2.895058e-04 1.072320e-04 1.141158e-03 
##          361          362          363          364          365 
## 4.935098e-04 4.322763e-03 4.292238e-03 5.731036e-04 1.015266e-03 
##          366          367          368          369          370 
## 1.222428e-02 4.520365e-04 8.836257e-04 4.301008e-06 1.025994e-03 
##          371          372          373          374          375 
## 1.155738e-03 7.250348e-05 6.368504e-04 2.767466e-03 2.247564e-03 
##          376          377          378          379          380 
## 1.409342e-04 8.629613e-05 3.288017e-06 1.437695e-05 1.175090e-06 
##          381          382          383          384          385 
## 2.532978e-03 7.736559e-05 4.522563e-03 9.199994e-04 1.636145e-09 
##          386          387          388          389          390 
## 9.179859e-03 1.939028e-03 1.048373e-03 2.006180e-03 9.043328e-05 
##          391          392          393          394          395 
## 2.043244e-04 1.819969e-04 5.282918e-03 1.047493e-05 3.865099e-04 
##          396          397          398          399          400 
## 1.595736e-04 9.715567e-07 1.649315e-03 4.660351e-04 2.702553e-06 
##          401          402          403          404          405 
## 3.265448e-04 2.975850e-05 1.782046e-04 6.869838e-03 1.578350e-03 
##          406          407          408          409          410 
## 1.320362e-03 6.239612e-05 3.516138e-03 6.107803e-03 5.268216e-03 
##          411          412          413          414          415 
## 8.229372e-06 1.461833e-02 9.274811e-04 1.289364e-05 2.072355e-05 
##          416          417          418          419          420 
## 2.455461e-05 3.639891e-04 2.562616e-04 2.381704e-04 5.211356e-04 
##          421          422 
## 3.704604e-04 1.143229e-03
#Create histogram
 
#A density plot of the residuals
plot(density(resid(model2))) 

#Create a QQ plot qPlot(model, main="QQ Plot") #qq plot for studentized resid 
leveragePlots(model2) # leverage plots

#Collinearity
vif(model2)
##   ydata$toptim  ydata$tmarlow ydata$tpstress 
##       1.279007       1.057358       1.330851
sqrt(vif(model2)) 
##   ydata$toptim  ydata$tmarlow ydata$tpstress 
##       1.130932       1.028279       1.153625
#Model 3 adding in gender 
#dummycode
ydata$gender=recode(ydata$sex,'0=1;1=2')
model3=lm(ydata$tpcoiss~ydata$toptim+ydata$tmarlow+ydata$tpstress+ydata$gender)
summary(model3)
## 
## Call:
## lm(formula = ydata$tpcoiss ~ ydata$toptim + ydata$tmarlow + ydata$tpstress + 
##     ydata$gender)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -34.785  -5.241   0.139   5.801  28.424 
## 
## Coefficients:
##                   Estimate Std. Error t value Pr(>|t|)    
## (Intercept)       56.17192    4.56967  12.292  < 2e-16 ***
## ydata$toptim       0.86391    0.11285   7.655 1.36e-13 ***
## ydata$tmarlow      1.03612    0.22525   4.600 5.62e-06 ***
## ydata$tpstress    -0.78439    0.08814  -8.899  < 2e-16 ***
## ydata$genderMALES  1.71591    0.90775   1.890   0.0594 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 8.991 on 417 degrees of freedom
## Multiple R-squared:  0.4449, Adjusted R-squared:  0.4396 
## F-statistic: 83.56 on 4 and 417 DF,  p-value: < 2.2e-16
plot(model3)

#Check assumptions
#Cooks distance
cooks.distance(model3)
##            1            2            3            4            5 
## 4.197679e-04 1.211857e-02 1.618401e-02 5.551629e-04 2.264980e-03 
##            6            7            8            9           10 
## 3.231304e-02 9.922411e-03 2.344161e-02 8.680495e-03 1.496076e-02 
##           11           12           13           14           15 
## 2.101761e-02 3.420758e-02 9.772502e-05 1.183614e-02 2.118482e-04 
##           16           17           18           19           20 
## 2.169598e-03 3.911430e-03 2.351318e-02 3.793327e-04 8.606866e-04 
##           21           22           23           24           25 
## 2.142289e-05 2.556939e-02 2.569283e-03 2.011169e-03 1.843008e-03 
##           26           27           28           29           30 
## 3.790730e-06 4.373871e-03 1.705208e-03 6.581347e-06 8.048003e-04 
##           31           32           33           34           35 
## 2.820298e-02 7.733624e-04 8.887682e-05 3.814116e-04 1.321204e-03 
##           36           37           38           39           40 
## 3.742400e-04 8.454544e-04 3.412412e-03 9.610690e-04 4.561408e-03 
##           41           42           43           44           45 
## 3.149815e-04 2.948345e-03 9.000356e-03 1.945060e-03 3.657561e-04 
##           46           47           48           49           50 
## 5.664595e-03 2.798421e-04 1.735660e-03 6.136839e-04 5.376345e-03 
##           51           52           53           54           55 
## 5.875351e-02 4.307659e-02 8.650980e-04 4.020425e-03 2.287152e-03 
##           56           57           58           59           60 
## 2.855970e-04 6.774712e-03 9.272567e-04 1.290323e-04 5.768160e-03 
##           61           62           63           64           65 
## 8.400061e-03 1.183437e-02 3.194868e-05 7.653155e-04 7.403000e-03 
##           66           67           68           69           70 
## 9.829533e-04 3.375666e-03 9.262198e-04 6.846940e-04 7.718024e-04 
##           71           72           73           74           75 
## 2.564665e-04 4.014454e-07 6.411572e-04 1.537949e-02 1.245155e-06 
##           76           77           78           79           80 
## 1.015045e-02 7.780475e-04 3.362903e-04 8.177173e-04 2.374995e-03 
##           81           82           83           84           85 
## 3.319909e-05 7.542612e-04 1.770498e-04 1.066118e-03 1.332752e-03 
##           86           87           88           89           90 
## 1.432210e-03 1.959148e-04 4.195844e-03 6.321568e-04 5.037242e-04 
##           91           92           93           94           95 
## 2.854249e-05 2.739247e-05 8.179369e-04 7.255416e-03 1.756023e-04 
##           96           97           98           99          100 
## 3.001939e-05 7.345742e-04 1.164994e-02 4.180816e-04 8.294681e-04 
##          101          102          103          104          105 
## 1.951387e-03 1.144373e-03 4.222578e-04 2.941887e-04 1.290234e-03 
##          106          107          108          109          110 
## 2.272432e-03 8.906600e-03 3.313974e-05 2.603343e-03 1.867603e-04 
##          111          112          113          114          115 
## 6.797506e-06 6.601685e-04 1.803555e-07 6.741639e-03 3.422203e-03 
##          116          117          118          119          120 
## 1.503357e-03 1.788088e-05 2.180131e-03 1.285808e-03 9.625718e-04 
##          121          122          123          124          125 
## 4.953154e-03 1.233853e-02 2.830108e-03 1.493332e-05 1.607792e-06 
##          126          127          128          129          130 
## 1.364789e-04 1.130918e-04 3.546110e-04 1.718545e-04 1.339119e-03 
##          131          132          133          134          135 
## 1.268833e-03 1.004576e-03 3.598073e-04 3.457427e-02 1.177440e-02 
##          136          137          138          139          140 
## 1.465423e-03 6.586219e-04 1.692150e-03 5.728757e-06 1.557841e-03 
##          141          142          143          144          145 
## 1.118464e-05 3.838866e-03 3.106583e-03 8.568775e-05 3.553826e-05 
##          146          147          148          149          150 
## 1.246649e-02 8.086788e-03 2.643602e-03 1.456764e-03 3.460046e-04 
##          151          152          153          154          155 
## 1.996944e-03 1.239686e-06 5.470631e-04 2.295693e-03 2.132351e-05 
##          156          157          158          159          160 
## 1.518723e-03 7.688412e-04 3.637548e-04 4.600586e-04 1.630315e-05 
##          161          162          163          164          165 
## 9.387885e-04 3.109898e-04 2.811556e-05 2.696406e-03 1.006388e-03 
##          166          167          168          169          170 
## 6.305138e-03 1.826157e-05 1.139847e-04 9.890679e-04 4.622274e-05 
##          171          172          173          174          175 
## 5.052888e-05 4.211214e-03 2.027140e-03 2.798842e-03 6.994851e-04 
##          176          177          178          179          180 
## 6.132901e-03 3.532210e-03 6.137388e-04 7.337603e-03 8.815952e-03 
##          181          182          183          184          185 
## 5.383469e-03 5.982091e-03 2.209876e-04 5.225183e-04 1.736916e-03 
##          186          187          188          189          190 
## 1.868019e-04 1.227234e-03 2.975273e-03 9.202072e-05 2.526248e-03 
##          191          192          193          194          195 
## 4.980183e-03 1.872972e-04 1.115488e-04 1.887004e-04 7.764134e-05 
##          196          197          198          199          200 
## 4.584208e-04 7.579627e-05 5.554306e-06 1.657696e-04 1.571302e-04 
##          201          202          203          204          205 
## 1.887963e-04 1.178353e-03 8.196134e-06 1.292070e-04 6.988422e-04 
##          206          207          208          209          210 
## 1.836346e-03 9.680570e-03 3.959528e-07 8.244072e-04 8.892333e-04 
##          211          212          213          214          215 
## 8.193632e-07 3.196681e-03 3.262213e-04 6.803670e-04 3.960848e-03 
##          216          217          218          219          220 
## 7.875614e-04 4.034885e-04 9.144393e-05 2.000980e-03 1.580416e-06 
##          221          222          223          224          225 
## 1.251735e-04 4.853038e-05 8.364939e-04 1.050282e-02 4.277987e-04 
##          226          227          228          229          230 
## 3.209268e-04 2.685046e-05 1.229495e-03 1.239696e-06 3.214409e-03 
##          231          232          233          234          235 
## 9.698840e-04 1.182820e-02 1.021546e-04 5.277873e-07 1.127162e-04 
##          236          237          238          239          240 
## 1.881968e-04 2.164420e-03 4.245212e-07 1.509754e-03 1.964100e-04 
##          241          242          243          244          245 
## 6.573517e-03 1.950828e-03 2.166570e-02 3.045035e-03 3.120855e-03 
##          246          247          248          249          250 
## 1.417851e-03 9.999129e-03 2.069329e-05 2.226214e-03 4.527933e-06 
##          251          252          253          254          255 
## 3.130802e-03 6.585050e-04 6.513967e-04 1.986348e-04 4.878517e-06 
##          256          257          258          259          260 
## 1.114737e-03 2.094699e-05 1.568792e-03 1.800409e-05 4.649115e-04 
##          261          262          263          264          265 
## 1.585284e-06 1.964130e-03 5.910648e-03 1.432358e-03 2.332845e-08 
##          266          267          268          269          270 
## 4.030685e-04 2.754178e-03 3.626445e-03 9.083886e-04 1.011470e-02 
##          271          272          273          274          275 
## 8.473261e-03 3.585344e-05 4.637712e-03 9.694950e-06 1.749809e-04 
##          276          277          278          279          280 
## 8.280766e-03 6.874586e-05 4.168752e-03 2.036204e-03 7.963922e-05 
##          281          282          283          284          285 
## 3.215642e-06 3.875906e-03 1.643751e-04 9.433776e-04 8.932770e-05 
##          286          287          288          289          290 
## 9.516064e-05 2.844443e-04 4.194760e-04 1.148010e-04 7.885892e-04 
##          291          292          293          294          295 
## 1.053573e-03 2.040621e-03 1.746390e-02 1.459593e-03 6.959472e-04 
##          296          297          298          299          300 
## 5.881585e-04 2.134193e-04 2.817285e-04 6.740716e-05 6.774610e-03 
##          301          302          303          304          305 
## 3.136208e-03 2.552473e-03 1.167534e-04 5.114699e-04 2.411268e-03 
##          306          307          308          309          310 
## 1.010636e-04 1.428834e-03 2.352240e-03 3.147816e-03 5.248036e-06 
##          311          312          313          314          315 
## 1.550405e-04 7.489822e-04 2.098045e-03 5.707228e-04 4.662413e-04 
##          316          317          318          319          320 
## 2.157255e-04 1.327133e-03 1.356481e-04 1.638573e-03 5.604195e-04 
##          321          322          323          324          325 
## 1.285783e-03 5.643174e-07 2.122820e-04 1.265711e-04 4.761796e-04 
##          326          327          328          329          330 
## 5.093441e-03 1.105760e-06 3.573770e-04 4.116176e-03 4.816687e-04 
##          331          332          333          334          335 
## 5.366328e-06 5.590812e-04 1.630920e-06 2.778076e-06 1.104570e-04 
##          336          337          338          339          340 
## 1.016314e-04 1.449896e-03 5.954211e-05 1.012711e-04 4.983936e-03 
##          341          342          343          344          345 
## 1.939200e-04 2.196996e-04 1.171223e-03 5.847426e-04 1.630944e-03 
##          346          347          348          349          350 
## 8.677706e-08 2.287167e-03 2.132634e-05 2.309821e-03 2.387297e-06 
##          351          352          353          354          355 
## 1.733529e-04 4.974783e-05 2.299280e-04 8.510232e-04 2.006994e-02 
##          356          357          358          359          360 
## 4.177580e-04 5.147529e-03 5.117771e-04 2.998917e-04 1.351073e-03 
##          361          362          363          364          365 
## 6.899045e-04 5.228499e-03 4.600990e-03 4.669527e-04 8.438978e-04 
##          366          367          368          369          370 
## 1.628354e-02 7.624859e-04 8.018057e-04 4.777453e-06 1.873800e-03 
##          371          372          373          374          375 
## 1.819473e-03 1.832969e-04 4.039127e-04 3.244130e-03 1.685125e-03 
##          376          377          378          379          380 
## 1.025939e-04 3.137450e-05 2.509090e-05 9.521722e-05 3.301360e-06 
##          381          382          383          384          385 
## 2.001649e-03 1.905834e-05 6.323092e-03 1.292753e-03 8.664254e-06 
##          386          387          388          389          390 
## 7.407007e-03 2.359656e-03 1.803314e-03 2.069562e-03 4.253547e-05 
##          391          392          393          394          395 
## 1.136323e-04 1.140534e-04 4.975520e-03 4.496540e-05 2.558134e-04 
##          396          397          398          399          400 
## 1.133500e-04 1.842323e-05 2.351419e-03 8.634116e-04 2.216123e-05 
##          401          402          403          404          405 
## 4.282417e-04 6.879309e-06 4.223486e-04 5.942391e-03 1.266087e-03 
##          406          407          408          409          410 
## 1.072475e-03 9.572363e-06 3.647850e-03 4.771934e-03 4.535670e-03 
##          411          412          413          414          415 
## 3.759026e-05 1.471750e-02 1.470953e-03 3.668430e-05 5.600332e-05 
##          416          417          418          419          420 
## 9.061783e-05 3.107488e-04 1.358236e-04 4.546927e-04 6.059712e-04 
##          421          422 
## 5.939650e-04 1.150430e-03
#Create histogram
 
#A density plot of the residuals
plot(density(resid(model3))) 

#Create a QQ plot qqPlot(model, main="QQ Plot") #qq plot for studentized resid 
leveragePlots(model3) # leverage plots

#Collinearity
vif(model3)
##   ydata$toptim  ydata$tmarlow ydata$tpstress   ydata$gender 
##       1.291821       1.075754       1.386133       1.050668
sqrt(vif(model3)) 
##   ydata$toptim  ydata$tmarlow ydata$tpstress   ydata$gender 
##       1.136583       1.037185       1.177342       1.025021
#Model 4 adding in interaction term gender*stress
#create interaction term
ydata$intgenstress=as.numeric(ydata$gender)*ydata$tpstress
model4=lm(ydata$tpcoiss~ydata$toptim+ydata$tmarlow+ydata$tpstress+ydata$gender+ydata$intgenstress)
summary(model4)
## 
## Call:
## lm(formula = ydata$tpcoiss ~ ydata$toptim + ydata$tmarlow + ydata$tpstress + 
##     ydata$gender + ydata$intgenstress)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -35.275  -5.186   0.235   5.836  28.198 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)        55.48477    4.86777  11.398  < 2e-16 ***
## ydata$toptim        0.86682    0.11318   7.659 1.33e-13 ***
## ydata$tmarlow       1.02637    0.22671   4.527 7.81e-06 ***
## ydata$tpstress     -0.69471    0.23450  -2.962  0.00323 ** 
## ydata$genderMALES   3.42931    4.24948   0.807  0.42013    
## ydata$intgenstress -0.06507    0.15766  -0.413  0.68000    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 9 on 416 degrees of freedom
## Multiple R-squared:  0.4451, Adjusted R-squared:  0.4385 
## F-statistic: 66.75 on 5 and 416 DF,  p-value: < 2.2e-16
plot(model4)

#Check assumptions
#Cooks distance
cooks.distance(model4)
##            1            2            3            4            5 
## 3.683116e-04 1.178676e-02 1.644426e-02 5.921484e-04 1.911109e-03 
##            6            7            8            9           10 
## 3.002889e-02 1.058362e-02 3.269252e-02 1.089707e-02 2.645437e-02 
##           11           12           13           14           15 
## 1.950181e-02 3.899282e-02 2.062466e-04 1.560389e-02 1.749112e-04 
##           16           17           18           19           20 
## 2.178498e-03 4.778911e-03 2.797210e-02 6.712949e-04 7.054063e-04 
##           21           22           23           24           25 
## 1.207878e-04 2.932034e-02 2.164039e-03 1.732762e-03 1.776822e-03 
##           26           27           28           29           30 
## 8.438375e-06 5.174328e-03 1.477976e-03 1.689095e-05 6.693014e-04 
##           31           32           33           34           35 
## 2.870488e-02 1.367956e-03 7.293801e-05 3.527067e-04 1.119045e-03 
##           36           37           38           39           40 
## 3.316852e-04 7.061671e-04 4.842655e-03 1.009167e-03 5.502894e-03 
##           41           42           43           44           45 
## 3.682338e-04 3.127571e-03 8.421178e-03 3.146750e-03 4.488304e-04 
##           46           47           48           49           50 
## 6.177542e-03 2.151253e-04 1.680959e-03 6.133201e-04 4.631203e-03 
##           51           52           53           54           55 
## 1.001518e-01 7.167624e-02 7.277119e-04 3.383246e-03 2.016119e-03 
##           56           57           58           59           60 
## 2.350115e-04 5.663261e-03 7.809555e-04 1.041914e-04 5.772658e-03 
##           61           62           63           64           65 
## 8.643961e-03 1.290567e-02 3.897609e-05 1.109233e-03 6.436888e-03 
##           66           67           68           69           70 
## 8.623124e-04 3.626219e-03 7.881128e-04 8.903274e-04 7.153656e-04 
##           71           72           73           74           75 
## 2.233232e-04 2.984345e-08 4.815933e-04 1.510601e-02 5.630290e-06 
##           76           77           78           79           80 
## 1.047456e-02 1.711374e-03 2.787357e-04 7.954411e-04 2.579801e-03 
##           81           82           83           84           85 
## 1.089499e-04 7.386730e-04 1.983170e-04 1.420885e-03 1.129245e-03 
##           86           87           88           89           90 
## 1.195295e-03 2.300407e-04 5.443813e-03 6.965125e-04 9.945682e-04 
##           91           92           93           94           95 
## 4.931185e-05 5.330733e-06 7.924720e-04 8.064176e-03 2.843402e-04 
##           96           97           98           99          100 
## 2.403436e-05 7.286398e-04 1.084999e-02 4.658209e-04 7.429898e-04 
##          101          102          103          104          105 
## 1.828286e-03 1.203475e-03 3.675654e-04 3.681127e-04 1.756676e-03 
##          106          107          108          109          110 
## 1.981103e-03 9.754931e-03 2.354276e-05 3.355610e-03 1.756504e-04 
##          111          112          113          114          115 
## 3.399336e-06 5.491710e-04 6.431870e-06 9.923681e-03 3.236055e-03 
##          116          117          118          119          120 
## 1.624381e-03 1.172468e-05 2.206231e-03 1.337213e-03 8.058774e-04 
##          121          122          123          124          125 
## 4.108480e-03 1.170423e-02 3.615455e-03 8.760133e-06 1.472871e-06 
##          126          127          128          129          130 
## 1.144466e-04 9.985317e-05 4.278511e-04 1.536422e-04 1.317295e-03 
##          131          132          133          134          135 
## 1.372503e-03 9.001797e-04 3.873307e-04 3.479622e-02 1.093127e-02 
##          136          137          138          139          140 
## 1.703842e-03 5.474786e-04 1.924143e-03 3.277167e-06 2.786034e-03 
##          141          142          143          144          145 
## 8.920424e-06 4.944636e-03 2.807581e-03 1.324993e-04 4.125021e-05 
##          146          147          148          149          150 
## 1.116912e-02 6.872314e-03 2.655669e-03 1.348741e-03 5.622060e-04 
##          151          152          153          154          155 
## 1.801882e-03 4.300273e-06 5.199456e-04 2.348412e-03 2.588558e-05 
##          156          157          158          159          160 
## 1.433510e-03 7.745445e-04 3.306319e-04 4.034340e-04 1.404673e-05 
##          161          162          163          164          165 
## 7.893235e-04 2.621461e-04 1.683682e-05 2.722357e-03 8.687618e-04 
##          166          167          168          169          170 
## 6.523981e-03 1.382915e-05 1.082414e-04 1.165182e-03 3.988461e-05 
##          171          172          173          174          175 
## 4.140182e-05 5.271212e-03 2.711268e-03 3.018363e-03 7.328011e-04 
##          176          177          178          179          180 
## 5.153066e-03 3.042363e-03 5.065314e-04 6.485450e-03 1.062737e-02 
##          181          182          183          184          185 
## 4.907789e-03 5.088367e-03 4.011871e-04 4.368638e-04 1.502905e-03 
##          186          187          188          189          190 
## 2.073765e-04 1.662333e-03 2.553593e-03 7.616832e-05 2.195589e-03 
##          191          192          193          194          195 
## 4.695773e-03 1.739894e-04 1.441452e-04 1.828951e-04 5.316878e-05 
##          196          197          198          199          200 
## 4.573369e-04 1.054271e-04 8.035693e-06 1.711428e-04 1.323312e-04 
##          201          202          203          204          205 
## 1.627571e-04 1.017564e-03 1.343025e-05 1.048791e-04 5.848871e-04 
##          206          207          208          209          210 
## 1.798831e-03 8.274625e-03 7.276398e-07 6.814074e-04 7.999788e-04 
##          211          212          213          214          215 
## 2.636695e-06 4.196595e-03 3.326340e-04 5.668468e-04 3.560563e-03 
##          216          217          218          219          220 
## 7.932861e-04 5.116409e-04 8.555473e-05 1.993805e-03 2.969881e-06 
##          221          222          223          224          225 
## 1.038617e-04 3.587097e-05 7.723367e-04 1.043787e-02 3.632631e-04 
##          226          227          228          229          230 
## 2.918906e-04 2.170598e-05 1.274902e-03 1.191147e-06 2.706881e-03 
##          231          232          233          234          235 
## 1.002980e-03 1.127884e-02 1.462014e-04 2.341247e-06 9.691864e-05 
##          236          237          238          239          240 
## 1.562080e-04 1.918112e-03 1.208516e-06 1.577807e-03 2.246451e-04 
##          241          242          243          244          245 
## 6.856845e-03 1.643265e-03 1.836275e-02 3.000762e-03 2.938401e-03 
##          246          247          248          249          250 
## 1.424394e-03 1.059559e-02 1.504406e-05 1.875266e-03 3.270413e-06 
##          251          252          253          254          255 
## 3.020612e-03 5.479769e-04 5.347724e-04 1.788675e-04 7.067421e-06 
##          256          257          258          259          260 
## 1.579329e-03 1.297818e-05 1.306728e-03 1.618743e-05 4.197649e-04 
##          261          262          263          264          265 
## 5.603899e-07 1.682458e-03 5.078914e-03 1.249098e-03 2.138688e-07 
##          266          267          268          269          270 
## 3.637804e-04 2.297152e-03 3.625006e-03 9.706161e-04 8.453090e-03 
##          271          272          273          274          275 
## 8.221500e-03 2.903182e-05 3.917443e-03 1.289357e-05 1.876943e-04 
##          276          277          278          279          280 
## 7.811538e-03 6.241586e-05 3.669516e-03 2.260326e-03 9.696378e-05 
##          281          282          283          284          285 
## 1.075501e-05 3.439622e-03 1.476867e-04 9.451931e-04 8.043664e-05 
##          286          287          288          289          290 
## 8.006605e-05 2.608237e-04 3.642120e-04 1.110590e-04 6.912347e-04 
##          291          292          293          294          295 
## 1.037926e-03 2.302189e-03 1.506446e-02 1.213464e-03 5.784579e-04 
##          296          297          298          299          300 
## 5.477090e-04 2.494791e-04 2.621553e-04 5.609695e-05 6.037795e-03 
##          301          302          303          304          305 
## 2.623062e-03 2.190796e-03 1.184799e-04 4.595411e-04 2.136557e-03 
##          306          307          308          309          310 
## 1.178762e-04 1.255140e-03 1.957706e-03 2.617487e-03 4.647854e-06 
##          311          312          313          314          315 
## 1.203079e-04 7.784063e-04 1.996550e-03 5.339854e-04 4.768896e-04 
##          316          317          318          319          320 
## 1.873916e-04 1.377090e-03 1.124370e-04 1.388773e-03 7.133237e-04 
##          321          322          323          324          325 
## 1.068561e-03 6.688153e-07 2.277190e-04 1.005425e-04 4.126710e-04 
##          326          327          328          329          330 
## 4.239403e-03 3.872964e-07 3.208166e-04 3.558723e-03 3.987121e-04 
##          331          332          333          334          335 
## 4.336547e-07 6.265775e-04 6.362877e-07 2.184722e-06 1.503205e-04 
##          336          337          338          339          340 
## 8.629488e-05 1.203576e-03 5.342945e-05 8.276439e-05 4.172393e-03 
##          341          342          343          344          345 
## 2.010043e-04 1.813412e-04 1.103757e-03 5.116216e-04 1.421534e-03 
##          346          347          348          349          350 
## 3.237523e-07 2.124335e-03 1.656522e-05 2.107369e-03 1.064613e-06 
##          351          352          353          354          355 
## 1.406517e-04 4.036877e-05 1.898873e-04 7.213794e-04 1.674341e-02 
##          356          357          358          359          360 
## 3.894949e-04 5.131893e-03 4.280909e-04 2.882945e-04 1.268574e-03 
##          361          362          363          364          365 
## 6.128398e-04 4.653312e-03 3.966617e-03 4.178548e-04 7.446793e-04 
##          366          367          368          369          370 
## 1.399947e-02 7.238569e-04 6.688734e-04 8.117083e-06 1.548044e-03 
##          371          372          373          374          375 
## 1.528814e-03 1.666879e-04 3.347214e-04 2.698099e-03 1.396226e-03 
##          376          377          378          379          380 
## 9.286380e-05 2.205523e-05 2.483743e-05 9.129923e-05 3.195493e-06 
##          381          382          383          384          385 
## 1.689474e-03 1.323054e-05 5.653313e-03 1.081766e-03 1.099372e-05 
##          386          387          388          389          390 
## 6.544350e-03 2.091592e-03 1.572124e-03 1.877792e-03 3.741160e-05 
##          391          392          393          394          395 
## 9.321803e-05 9.272344e-05 4.160715e-03 4.927077e-05 2.306507e-04 
##          396          397          398          399          400 
## 9.357719e-05 2.022122e-05 1.968139e-03 7.278751e-04 2.329951e-05 
##          401          402          403          404          405 
## 3.854759e-04 5.634831e-06 3.591673e-04 4.942919e-03 1.053912e-03 
##          406          407          408          409          410 
## 9.012352e-04 9.246236e-06 3.098641e-03 3.988597e-03 3.768479e-03 
##          411          412          413          414          415 
## 3.350123e-05 1.248138e-02 1.320237e-03 3.436121e-05 4.610486e-05 
##          416          417          418          419          420 
## 7.284765e-05 2.678897e-04 1.200142e-04 3.945567e-04 5.042197e-04 
##          421          422 
## 5.055682e-04 9.536790e-04
#Create histogram
 
#A density plot of the residuals
plot(density(resid(model4))) 

#Create a QQ plot qqPlot(model, main="QQ Plot") #qq plot for studentized resid 
leveragePlots(model4) # leverage plots

#Collinearity
vif(model4)
##       ydata$toptim      ydata$tmarlow     ydata$tpstress 
##           1.296847           1.087557           9.791870 
##       ydata$gender ydata$intgenstress 
##          22.979404          27.390333
sqrt(vif(model4)) 
##       ydata$toptim      ydata$tmarlow     ydata$tpstress 
##           1.138792           1.042860           3.129196 
##       ydata$gender ydata$intgenstress 
##           4.793684           5.233577